DIFFERENTIAL SYSTEMS WITH FUCHSIAN LINEAR 
PART: CORRECTION AND LINEARIZATION, 
NORMAL FORMS AND MULTIPLE ORTHOGONAL 

POLYNOMIALS 



Abstract. Differential systems with a Fuchsian linear part are 
studied in regions including all the singularities in the complex 
plane of these equations. Such systems are not necessarily ana- 
lytically equivalent to their linear part (they are not linearizable) 
and obstructions are found as a unique nonlinear correction after 
which the system becomes formally linearizable. 
More generally, normal forms are found. 

The corrections and the normal forms are found constructively. 
Expansions in multiple orthogonal polynomials and their gener- 
alization to matrix-valued polynomials are instrumental to these 
constructions. 



1.1. The problem. A first order hnear system 



is caUed Fuchsian if all its singularities in C are regular. For such 
systems A{x) has the form 



and then ([T]) is singular ai po, . . . ,Ps+i, oo, all of them regular singular 
pointsEH This linear context appears in a wide range of problems of 
mathematics and physics, and has been the topic of extensive studies. 
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1. Introduction 



(1) 




(2) 




See H5.1l5.2l for a classification of singularities. 

^The peculiar choice of the limits of summation in ^ is for the convenience of 
the present paper. 
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The present paper studies nonlinear systems having a Fuchsian hnear 
part: 

(3) ^ = A{x)u + -^{{x,u), ueC' 
ax Q{x) 

with A{x) given by ([2]). The term collects the nonlinear terms in u 
and Q denotes the product 

5+1 

(4) Q{x) = ll{x-p,) 

j=0 

The question addressed is: under which conditions systems ([3]) are 
analytically equivalent to their linear part ([1]) on a domain in x large 
enough to encompass all the singularities po, . . ■ ,Ps+i? 

This is a question of classification of differential equations in singular 
regions. It is well known that in a neighborhood of a regular point, 
all equations are analytic equivalent. But this is no longer true on a 
domain including singularities. 

1.2. Motivation. The question of linearization and, more generally, 
of equivalence, is a fundamental problem in the theory of differential 
equations. 

Besides the clear theoretical interest, linearization and equivalence 
are used as instrumental methods in control theory, and in devising 
algorithms for numerical and symbolic calculations (see, to cite just a 
few authors, [12], [13], [E], [IZ], [IB])- Often the methods used are 
developed from the method of equivalence introduced by Cartan [2] 
to decide whether two differential structures can be mapped one to 
another by a transformation taken in a given pseudogroup [20], [25] . 
In the case of differential equations, the method was used for regular 
systems and extended for a neighborhood of one singular manifold [21], 
[22]. 

There are many classes of problems that reduce to equations ([2]), 
([3]). Vector fields with an eigenvalue equal to 1 at a singular point 
can be reduced to the study of a Fuchsian system near one singularity 
by eliminating time. Similarly, vector fields in regions containing two, 
or more singular points are reducible to Q. More generally, Hamil- 
tonian systems with polynomial potentials near particular, periodic, or 
doubly-periodic solutions can be reduced to ([3]) [1], [5]. 

In dimension one nonlinear equations near a singular point (regu- 
lar, or irregular) were studied in detail by Martinet and Ramis. Au- 
tonomous nonlinear equations near a singular point have been also 



DIFFERENTIAL SYSTEMS WITH FUCHSIAN LINEAR PART 



3 



studied, have generated deep results, and are relatively well under- 
stood. 

The present paper considers systems in regions encompassing two, or 
more, singularities. The study also reveals a close connection between 
three concepts: linearizability, integrability and multiple orthogonal 
polynomials, and leads to generalizations of the latter. 

The result of Theorem[3] has an interesting similarity with needed 
corrections found in other problems: Ecalle and Vallet showed that 
resonant systems are linearizable after appropriate correction ; also 
Gallavotti showed that there exists appropriate corrections of Hamil- 
tonian systems so that the new system is integrable [I9], convergence 
being proved later by Eliasson [I6j. This suggests the possible existence 
of an underlying general structure. 



1.3. Notations used. 

Denote N= {0,1,2,...}. 

For u, A G C"' and m G N'^, denote 

|m| = nil + ^2 + . . . + rrifi, 
m ■ A = miAi + . . . + mdXd 

A^d(C) denotes the d x d matrices with complex entries. 
C^[a;] denotes the polynomials in x with coefficients in C^. 
A^TvN denotes the polynomials in x with coefficients in Ai^i^C). 
Qi denotes the unit vector in with coordinate i equal to 1. 



1.4. Prior results. 



1.4.1. Systems near one singularity. If an equation ©-Q is studied 
for X in a domain small enough to include only one of the singularities 
Pj, and placing this singularity conventionally at a; = 0, such a system 
has the form 

(5) '^ = -Lix)M + -Ux,u), ugC' 

dx X X 

with L{x) a matrix analytic at and f (x, u) analytic at x = 0, with a 
zero of order two at u = 0. 

If the eigenvalues of L{Q) are linearly nonresonant, in the sense of 
(169|) . then after an analytic change of coordinates the matrix L{x) can 
be assumed constant L{x) = L(0) = L (see §5.2p . 

Theorem[T] shows that analytic linearization near one singularity holds 
for generic systems [6]: 
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Theorem 1. Consider an equation 




Lu -\ — f (x, u) 



X 



where f (x, u) is analytic for \u\ < R (R > 0) and for x E D, where D 
is a disk centered at (or an annulus centered at Q). 

It is assumed that L is a diagonalizahle matrix with eigenvalues 
Ai, \d satisfying the following diophantine condition: there exist con- 
stants C,!^ > so that 



for all i G {!,...,(/}, m G N'^, |m| > 2, and for A; G N (respectively, 
k G Z). 



Then the system ^ is analytically equivalent to its linear part 



for X e D' d D and |w| < R' < R. 

Furhermore, D' can be made arbitrarily close to D if R' is small 
enough. 

Also, the analytic equivalence map is unique if no Aj is an integer. 

As a clear consequence of Theorem[T] the study of the local analytic 
properties of systems ([HD reduces to the study of the linear ones 

One application is to the study of integrability. Linear equations 
(IH]) have a full number of linearly independent first integrals which are 
analytic for small x and w, except for the singular manifolds x = 
and/or w = 0, where (generically) these first integrals have accumula- 
tion of poles: they are not meromorphic [6]. Then the same holds for 
the nonlinear systems ([6]). 

A posteriori, it is clear that one should expect that first integrals, 
when they exist, are usually not meromorphic. On the other hand, 
when looking for first integrals in concrete problems, meromorphicity 
is often assumed. Such assumptions simplify the study, since under 
(some form of) regularity assumptions the integrals are guaranteed to 
have specific types of expansions (e.g. Poincare series, Laurent series) 
allowing for a study by successive approximations. 

1.4.2. Systems in regions with two regular singularities. As seen in 
§1.4.11 near one singularity generic systems ([2])-([3]) are linearizable, and 
integrable. Let us consider a larger domain of x, so that it includes two 



(7) 



A; + m-A-A, > C (|m| + |A;|) 



V 




dw 

dx 



—L w 

X 
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singularities. 



More precisely, consider ([3]) for 5 = 0: 



dvi f 




f (x, u) 



Of course, the position of the two singularities po, pi can be arbitrarily 
placed in the complex plane using a linear change of variables. 

The situation changes dramatically from the study near one singular- 
ity ([6]): equations (Q are not necessarily linearizable for x in a domain 
D containing both singularities po and pi. 

This fact can be understood by applying Theorem[l]near singularity 
Po, and then, near singularity pi. Then (genericallju) there exists a 
unique linearization map analytic ai x = Pq, and there exists a unique 
linearization map analytic at a; = pi. But the two maps do not nec- 
essarily coincide, as the map analytic at x = pj is usually branched at 
x = Pi-j {j = 0,1). 

This fact has important consequences, as it turns out that lineariz- 
ability and integrability are intimately connected [7] : 

Theorem 2. Consider an equation ^ in dimension 1: 



If equation [W\] is not analytically linearizable then no single-valued 
integrals exists on domains encircling ±1 for generi^ aQ,ai. 

It is then important to determine which equations are linearizable, 
and which are not, and to find the equivalence classes with respect to 
analytic equivalence in this semi- local context. 

This problem was studied in [5] for S = 0; the results follow from 
Theorem[3] proved here. 



2.1. Setting. Let S* > 0, let po,pi, . . . ,ps+i be distinct points in C 
and let Aq, Ai, . . . , As+i he d x d matrices. 

Consider the differential system ([3]) with ([2]), (jlj). 

Let D be a simply connected domain in the complex plane containing 
the singularities Po, . . . ,Ps+i- 

"'Existence of an analytic lienarization map for x close to pj follows from Theo- 
rem[T]if Aj is diagonalizable, satisfies the Diophantine condition ([7]), and is linearly 
nonresonant, ^5.3\ 

^Precise conditions are given in [7j. 



(10) 




2. Main results 
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It is assumed that f (x, ■) has a zero of order two at u = 0, and is 
assumed analytic for x & D and small |u|. 

The denominator Q of the nonlinear part of ([3]) simply shows that 
this nonlinear part may have at most poles of order one at Po? • • • iPs+i- 

Denote 



5+1 

'^3 



(11) Aoo = Y. ^3 

Note that A[x) ~ A^x~^ as x ^ oo. 

Definition: A matrix M G A^rf(C) is called nonlinearly nonresonant 
if its eigenvalues Ai, . . . , satisfy 

(12) A; + A ■ m - Ai ^ 0, for all A; G N, m G N"', |m| > 2, z = 1, . . . , ci 

Assumption: It is assumed that each matrix Aq, . . . , v45+i, 
is nonlinearly nonresonant. 

2.2. Main results. 

Theorem 3. Consider an equation ^-^) under the notations and 
assumptions of 

Then there exists a unique correction of the nonlinear part as a for- 
mal series 

(13) 0(X,U)= Yl ^mla^)"'" 

meN<^, |m|>2 

where 0ni(^) ^^^^ polynomials in x of degree at most S, such that the 
corrected system 

(14) ^ = A(x) u + [f (x, u) - 0(x, u)] , ugC^ 
ax Q[x) 

is formally linearizable through a change of variable as a series 

(15) u = w + h(x, w) = w + hm(x)w™ 

|m|>2 

with hin(x) analytic on D. 
Remarks. 

1. Equation ([3]) is formally linearizable if and only if (f>{x,u) = 0. 
Therefore this unique correction gathers the obstructions to lineariz- 
ability. 

2. Obviously, if an equation is not formally linearizable, then it is 
not analytically linearizable either. 
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Since equations ([3]) are not necessarily linearizable, then they are 
not all equivalent either. TheoremU] provides the classification of these 
equations by specifying formal normal forms: 

Theorem 4. Consider an equation under the notations and 

assumptions of ^2.1\ 

Then there exists a unique formal series 

(16) V'(x,u)= ^m{x)n"' 

meN'*, |m|>2 

where ip^{x) are polynomials in x of degree at most S, such that the 
system ^ is formally equivalent to 

(17) ^ = 4Ww + ^V(x.w), wee' 

through a change of variable as a series of the form ( 173]) with hm{x) 
analytic on D. 

Comments about convergence. While 0j„(x), V'm(^) PolY" 
nomials, and hin(a;) are analytic functions, the series in w f|T3l) . f|T5|) . 
f[TB]l are formal. 

Convergence of f[T^ and of the linearization f[T^ were proved in 
[9J for 5 = under the supplementary assumptions that Aq, Ai are 
simultaneously diagonalizable, and have eigenvalues with positive real 
parts. 

For convergence to hold in general it is reasonable to expect that the 
nonresonance condition (fT2l) should be strengthened to a diophantine 
condition. This could be a Siegel type condition, like ([7]), or even a 
Bruno type condition; the latter was proved by Yoccoz to be optimal 
for the problem of iteration of holomorphic germs in one variable, and 
under this condition convergence holds in the linearization problem for 
vector fields at a resonant point, as proved by Ecalle and Vallet [T5] . 



3. Proofs 

3.1. The structure of the proofs. The main steps in the proofs of 
Theorems [3] andH] follows. Using power series expansions in 

w it turns out that the terms hm(x) in (fT5|) satisfy recursive systems 
of differential equations which are first order, linear nonhomogeneous. 
Remarkably, their linear part is Fuchsian. The burden of proof is car- 
ried by the study of these Fuchsian nonhomogeneous systems, and the 
structure of their solutions turns out to be rich and interesting. 
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3.2. A fundamental lemma. Consider a Fuchsian equation with a 
nonhomogeneous term: 

(18) y'(x)+ B{x) y{x) = -i- g(x), y, g G C^, B{x) G Mn 

Q[x) 

where 

s+i , s+i 



(19) S(a;) = ^-__i?^, Q{x) = \{{x~p,) 

Let 



a; — Pi 

3=0 j=0 



5+1 



(20) B^ = J2Bj 
It is assumed that 

(21) k + Bj are invertible for all A; G N and j = 0, 1, . . . , S" + 1, oo. 

Lemma 5. Consider the equation l[T^) - l[T^) under the assumption l[21\) . 

Let D be a simply connected domain containing the singularities 
Po, . . . ,P5+1- 

Then for any function g(x) analytic on D there exists a unique 
0(x) G C^[a;] polynomial of degree at most S so that the corrected 
equation 

(22) y'{x) + B{x) y{x) = ^ [g{x) - 0(a:)] 

has a solution y{x) which is analytic on D. 

The proof of Lemma[5]is the topic of ^ Before providing its details, 
which involve both analytic and algebraic arguments, and leads to gen- 
eralizations of multiple orthogonal polynomials to matrix-valued ones, 
let us assume Lemma[5]true for the moment and show its usefulness by 
providing the proofs of Theorems^] andlH 

3.3. Proof of TheoremlSl The change of variables (JTSll provides a 
linearization of (1141) if and only if h satisfies the nonlinear partial dif- 
ferential equation 

(23) d^h + d^h Aw = Ah+ -i- [f (x, w + h) - 0(x, w + h)] 

Q{x) 

Searching for solutions of fl2^ as power series in w, denote by h„ the 
homogeneous part degree n of h(x, w); in the notation ([T^ we have 

(24) h„(x,w)= hm(x)w'" 

|m|=n 
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We use similar notations to denote by and (j)^ the homogeneous part 
of degree n of f , respectively 0. 

Equation (!23|) splits into blocks of systems of ordinary differential 
equations for {hm}|m|=n, recursive in n: 

(25) 9^h„ + dwhnv4w-ylh„ = — ^R„(x,w), n>2 

Q{x) 

where 

(26) R„ = f„ - 0„ + R„ 

with R„ a polynomial in <p^, h^, with |m| < n, and R2 = 0. 

Each h„ and 0„ are to be determined from (!25|) - (!26|) inductively on 
n. This can be done as follows. 

Note that, remarkably, for each n, fl25p is a Fuchsian nonhomoge- 
neous system. 

To see this, denote by Vn the space of vector- valued polynomials in 
w G C^, homogeneous degree n: 

(27) P„ = I q ; q(w) = J] q^w"", q^ G 

meN'',|m|=n. 

Note the canonical basis of the linear space Vn'- 

(28) rin,j = w'^e,, |m| = n, i = 1, . . . , d 

and its finite dimension dimVn = N = N{d,n)^ 

Denote by B{x) the linear operator on Vn taking to 

B{x)hn = dwh„ Aw - Ahn 

We have 

5+1 ^ 

(29) B(x) = B. where Rq = dwq^jW - A,:q 

^-^ X — Pi 

j=0 

hence in the canonical base fl2Sl) B{x) is given by a Fuchsian matrix: 
for each n the system (!25|) . (126|) has the form (l22l) (with g = f„ + R„ 
and (f> = (f)n)- 

Note that -Boo = dwq^ooW — Aooq. 

Theorem[3] now follows from Lemma[5]once we show that its assump- 
tion (I2T!) is satisfied. This is proved by the following result: 



^In fact N ^d{n + d- l)!/n!/(d - 1)!. 
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Lemma 6. Let M be a d x d matrix with eigenvalues Ai, . . . , (not 
necessarily distinct). 

Let Jm he the following linear operator on 

(30) ( J^q) (w) = (dwq) Mw - Mq 

Then the eigenvalues of Jm are \ ■ m — Xi, i = 1, . . . ,d, |m| = n. 

The proof of Lemma[6] is found in §5.41 

In the notation fl29l) . fl30|) we have Bj = Jaj for all j = 0, . . . , 5+1, oo. 
Therefore, in view of LemmaEl (fT2l) implies (12T|) . □ 

3.4. Proof of Theorem[4l The change of variables (fT5l) transforms 
(IT^ into fll7p if and only if h satisfies the nonlinear partial differential 
equation 
(31) 

d^h + d^h ^Aw + ^V'(a;, w)^ = Ah + ^ [f{x, w + h) - ?/?(x, w)] 

Expanding in power series in w we obtain, as in the proof of Theo- 
rem[3l that the homogeneous parts h„ (n > 2) of h satisfy systems of 
equations which can be solved recursively in n: 

(32) 9^h„ + dwh„Aw-ylh„ = — ^R„(x, w), n>2 

Q{x) 

where 

(33) R„ = f„ - V„ + Rn 

with R„ a polynomial in ij)^, hm, dwhm, fm with |m| < n, and R2 = 0. 

Each and can now be determined from ( !32l) as in the proof of 
Theorem[31 □ 

4. Proof of Lemma[5] 

The proof consists of several parts. Section §4.21 shows that if the 
correction (f) exists, then it is unique. 

Section §4.31 proves LemmaOfor the case when g(x) is a polynomial. 
The proof is algebraic, and solutions are found as expansions in terms of 
multiple orthogonal polynomials and their generalizations introduced 
here. 

Section §4.41 proves Lemma[5] for the case when the eigenvalues of all 
matrices Bj have positive real parts. An analytic approach is used. 

Section §4.51 shows that the general case can be reduced to §4.4[ thus 
completing the proof of LemmaO The algebraic results of §4.31 are 
used to construct this reduction. 
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4.1. Existence of local analytic solutions. The existence and unique- 
ness of an analytic at pj solution of f[T5]) follows from the rest of ^ 
Due to its own intrinsic interest the statement and a simple proof are 
included here. 

Lemma 7. Under the assumption l[21\) equation l[T^) - l[T^) has a unique 
solution y(x) which is analytic at the singular point x = pj (for < 
J<S + 1). 

Proof of Lemma\7^ 

Substituting a power series y(x) = yo + (x — pj)yi + (a; — Pj)^y2 + . . . 
in (ITS]) one obtains a recursive system for the coefficients yo, Yi, y2, • • 

BjYo = go 

(1 + Bj)y^ 

{k + Bj)yk = rfe(yo, . . . , y^-i, g) 

which has a unique solution by the assumption fl2T|) . 

The series converges since equation f|T8|) is linear nonhomogeneous. 

□ 

4.2. Uniqueness of corrections. If there are two corrections (f> = <pi, 
and (p = </>2, so that equation fl22l) has the analytic on D solutions 
y = yi, respectively y = y2, then y = yi — y2 is an analytic solution 
of f|T8|) with g = 02 ~ 01 3^ polynomial degree at most S. Lemma[9] 
shows that this implies (pi = 02 ^.nd yi = y2. 

The following result is needed. 

Lemma 8. If g{x) is a polynomial then any solution of (T^-^M) is 
polynomially bounded for x oo. 

Proof of Lemma^ 

Denote by W{x) the integrating factor of f[T5]) : an invertible matrix 
satisfying 

(34) ^Wix) = W(x)B(x) 

ax 

Note that W{x) is a fundamental matrix for the equation 

(35) v = v5(x) 

and that W is the inverse of a fundamental matrix of the linear equation 
y' + 5(x)y = 0. 
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Since for ( l35l) the point at oo is Fuchsian the matrix W can be 
chosen such that W{x) = where B^o is a constant matrix 

whose eigenvalues (not counting the muhiphcity) are also eigenvalues 
of Boo, with $(a;) analytic at oo and $(x) = / + o(l) as x ^ oo (see 
§5.21 for details). Therefore 

(36) W{x) = x^°^ (J + o(l)) for X ^ oo 
The general solution of (fTSj) is 

(37) W{x)-^$, + W{x)-^ j Q{t)-'W{t)g{t)dt 

where the path of integration avoids the singularities pj, the integral 
starts at an arbitrary, but fixed point, and ^ is a constant vector. 
Denote g(a;) = gmX"' + . . . + gix + go. 

For X — > oo we have W{x) ~ x^°°, Q{x) ~ x'^"'"^ and g(x) ~ gmx"^ 
therefore ( !37l) is polynomially bounded for x — oo. □ 

Lemma 9. Consider equations [IE)-[T^) under the assumption (21\) . 

If g(x) is a polynomial of degree at most S then equation [T^) has 
a solution which is analytic at all the singular points po, . . . ,Ps+i only 
for g(x) = 0, and in this case the only analytic solution is y(x) = 0. 

Proof of Lemmo^ 

Assume that there exists y(x) a solution of (1181) which is analytic 
at all the singular points po, . . . ,Ps+i of the equation. Since (ITSl) is 
a linear nonhomogeneous ordinary differential equation then y(x) can 
have no other singularities and therefore y(x) is entire. By Lemma[8] 
then y(x) is a polynomial. 

If y(x) = then f|T8|) implies g(x) = 0. 

Assume, to obtain a contradiction, that y(x) ^ 0. Denote y(x) = 
y^x'' + . . . + yo with y„ ^ and g(x) = g^x'" + . . . + gix + go. 
Substituting these expansions in ([T5]l . for x — * oo the left-hand-side 
has the dominant term {n + B^)ynX'^~^ and the right hand-side has 
the dominant term g^x'""'^"^. Since m < S then m — S — 2<n — 1 
therefore we must have (n + -Boo)yn = 0, and by assumption fl2T]) then 
y„=0, which is a contradiction. □ 

4.3. Algebraic structure: proof of Lemma[5] for g(x) polynomi- 
als using expansions in multiple orthogonal polynomials and 
their generalizations. 

In this section it is assumed that g(x) are polynomials. LemmaO 
stating existence of corrections and of corresponding analytic solutions 
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(which turn out to be polynomials in this case) is proved using ex- 
pansions in multiple orthogonal polynomials and their generalizations, 
introduced here. 

Let B{x) be a Fuchsian matrix (fT9|) . with the notation (!20|) . Condi- 
tion is not used until UX2\ 

Let W{x) be a fundamental matrix of (!35|) . 
For any n = 0, 1, 2, . . . denote 

(38) n = {S + l)m + i with m = [n/{S + 1)J , and z = 0, 1, . . . , 5 

A precise notation requires m = m„ and i = in, but the subscript 
will be omitted to avoid a heavy notation. Anyhow, the one-to-one 
correspondence n (m, i) implied by is only used in §4.3[ 
With the notation fl38|) . define 

(39) P!,'^Hx)^Pn{x)=W{xr'^ [x^Q{xrW{x)] 
Remarks. 

1. In dimension one (A^ = 1) we have Bj = bj G C, and W{x) has 
the form 

5+1 

(40) W{x) = Yl [X — Pj) ^ 

j=0 

la. In particular, for S = only i = can occur in ( l38l) . and ( !39l) is 
the Rodrigues formula defining the Jacobi polynomials Pn°'^^\x) (up 
to a multiplicative constant). 

lb. For S = 1 formula (155]) defines the Jacobi-Angelesco multiple 
orthogonal polynomials introduced and studied in [21] (see also [23] )• 

Ic. For higher 5* the Pn{x) defined by fl39l) are a natural extension 
of the Jacobi-Angelesco polynomials. 

2a. For higher dimensions (A^ > 1) and 5 = formula fl39|) was stud- 
ied in [1^ where it is shown that Pn{x) are polynomials and that they 
form a complete system. In addition they satisfy a two step recurrence 
relation, and, in the commutative case, other beautiful properties of 
orthogonal polynomials also hold: they are eigenfunctions for a second 
order differential operator, and they are orthogonal with respect to the 
weight W{x) (when integrals exist). 

2b. For 5 = the polynomials occurring in the proof of Theorem[3] 
were studied in [11]. 

For the proofs in the present paper we need to show that formula (1391) 
defines polynomials indeed (done in Proposition fTOl) . that they form a 
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complete system (Propositiondl]), and that ( 122|) is solved in terms of 
them, 

4.3.1. Pn{x) are polynomials. Note first that Pi{x) = for z = 0, 1, ... 5*. 

The following reformulation of Rodrigues-type formula (J39ll helps 
simplify the proofs considerably (the case S = appears in [ID]). 

We let ^ = d^. 

ax 

Proposition 10. Denote by Ak the following linear operators onM.N[x]: 

(41) Ak = kQ'{x) +Q{x)B{x) + Q{x)d, 
Then the functions Pn defined by (3^) satisfy 

(42) Pn{x)=AiA2...Amx' 
for n > S + 1. 

As a consequence, Pn{x) are polynomials in x (matrix-valued). 
Proof of PropositionlTU 

Calculating successively the m derivatives in (13^ we obtain, using 
(El and dUD, 

(43) ^ Ix'Q'^W] = mx'Q'^-^Q'W + x'Q'^W +Q"'Wd^x' 
dx '- -' 

= Q'^-^W [[niQ' + QW-^W) x' + Qd^x'] = Q'^-^WAmx' 

then 

(44) ^,WQ-W] =^[Q-'WAr.x^] 

= Q^-^WAn,.iA„^x' 
and so on. □ 

Proposition 11. (i) The polynomial Pn defined by [SU) has degree at 
most n, and the coefficient of x" is 

m 

(45) Cn = X[{j+n + B^) forn>S + l 

i=i 

(a) If Boo satisfies 

(46) k + i?oo is invertibe for all k > S + 2, /c G N 

then ( f73] j is an invertible matrix, and therefore the degree of Pn is 
exactly n. 
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As a consequence, the set {P„(x)}„,gN is complete in the sense that 
any polynomial i{x) G C^[x] can be written as an expansion in Pn's: 

degf 

(47) f (x) = J2 Pn{x%, for some f„ G 

n=0 

and the expansion is unique. 
Proof of Propositionni\ 

(i) Retain only the dominant coefficients in fj^Tl) using the fact that 
Q{x) ~ x^^"^ and B{x) ~ x'^B^o as a; — > oo. We then have 

Ak ~ x^+i \k{S + 2) + B^ + 

Therefore 

Amx' ~ [m{S + 2) + B^ + i] 

then 

[{m - 1){S + 2) + B^ + I + S + I] [m{S + 2) + B^ + i] x'+^^^+^^ 
and so on. Using (1421) therefore 

m 

Pn{x) r^\[[j{S + 2) + B^ + i + {m- j){S + 1)] 
i=i 

which equals C„x", proving 

(a) Condition (1461) clearly ensures that fHSl) is an invertible matrix. 
Invertibility of the dominant coefficient of P„ implies, using induction 
over the degree of f, the existence of the unique expansion (jlT]). □ 

4.3.2. Proof of LemmalM for g{x) polynomial. Assume fl2Tl) satisfied. 
Denote W = ^W. Note that W is a fundamental matrix of solutions 

for the Fuchsian equation v' = vB where B = B — Q'/Q. Note that 
B = —^B~ with Bi = Bi — J, and therefore 



(48) Soo = Soo - (^ + 2)1 

Consider the polynomials associated to by a formula (1391) : P, 
P„. Note that we have 

(49) P^^s+i{x)=Q{x)W{x)-'^-^ [x^Q{xrW{x)] 



(W) 
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Using the integrating factor W equation ( l22l) can be rewritten as 

(50) ^ [ W{x)y{x) ] = W{x) [g(x) - (/)(a;)] 

In view of ([311), (Si) and ([SnD we see that if g = Pn+s+i then y = P„ 
is a solution of (|T8l) . 

For general polynomial g we only need to expand g in the special 

polynomials -P„, as guaranteed by Proposition[TT] (note that fl2Tl) . PHI) 

imply that Boo satisfies (H6l) ): 

deg g 

g(^) = Yl ^"'^^^ 
n.=0 

Letting 

deg g s ^ 

y{x) = X] ^«-5-i(a;)gn and 0(x) = ^ P„(a;) g„ 

n.=S'+l n=0 

the polynomial y{x) satisfies fl22|) . □ 
Remark. Note that degy(a;) =degg(a;) — S — 1. 

4.4. Analytic structure: proof of Lemma[5] when the eigenval- 
ues of all matrices Bj have positive real parts. In §4. li the unique 
solution of ( fTSj) which is analytic at x = po was constructed via its Tay- 
lor series. LemmafT2] shows that if the eigenvalues of all matrices Bj 
have positive real parts this solution can be written in an integral form 

Under the assumptions of §4.41 we can choose the starting point of 
the integral in (137]) to be a singularity of the equation, and the general 
solution of (fTSj) has the form 

(51) W{xy^^ + y{x) 
where 

(52) y{x) = W{x)-^ I Q{tY^W{t)^{t)dt 
and i e C^. 



Lemma 12. If the eigenvalues of Bq have positive real parts, then f5E 
is the unique solution of [T^) which is analytic at x = Pq. 

Proof. 

Let us first convince ourselves that (l52l) is analytic at a; = po- 
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There exists a fundamental matrix for ( 135|) of the form ^^0(2;) = 
{x — po)^"^o{x) where $0(2^) is analytic at po, $o(Po) = -^5 and i?o is a 
constant matrix (see §5.21 for details). Therefore W{x) = TqWo{x) for 
some constant, invertible matrix Tq. Note that (152|) is the same if we 
replace W{x) by Wo{x). 

Expanding in a convergent series, 



00 



Q(t)-%(t)g(t) = )_^(t-po)"fn 

n=-l 

we see that it is enough to show that (x — Po)~^° f^^ (t — po)^°~^"' dt is 
analytic at Pq for all n > —1. This is clearly true: the usual integration 
formula for powers applies because the matrix Bo+n+l is invertible due 
to the assumption (12T1) (and any eigenvalue of Bq is also an eigenvalue 
of Bo). 

Uniqueness of the analytic at po solution of 0181) was proved in §4.11 
Alternatively, a direct proof is by noting that W{x)~^$, is analytic at 
Po only for the constant vector ^ = 0. Indeed, assuming otherwise, 
W{x)^^^ is a nonzero solution of y' + B{x)y = which is analytic 
at X = Po and this implies that —Bq has an eigenvalue in N, which 
contradicts assumption fl2T]) . □ 

Note that in Lemma[T2]we can replace po, Bq by any pj, Bj, and these 
results are now used to determine the condition for the solution fl52l) 
to be analytic at other points pj as well. Rewriting fl52l) as 

(53) 

y{x) = W{x)-^ Q{t)-'W{t)g{t)dt+W{x)-^ Q{t)-^W{t)g{t)dt 

and using Lemma[T2] it follows that y(a;) is analytic at the point pj if 
only if 

(54) r Q{t)-^W{t)g{t)dt = 

•Jpo 

Then y(x) is analytic on D if flS^ holds for all j = 1, 2, . . . , S* + 1. 
Thus a correction (f){x) needs to be determined so that 

(55) r Qity^Wit) [git) - cj>{t)] dt = 0, Vj = 1, 2, . . . , ^ + 1 

•Jpo 



Denote <p{x) = Yld=o4>i^^- Conditions (155!) become 



5 

^/ x'Q{x)-^W{x)dx4)i= Q{x)-^W{x)g{x)dx, J = 1,...,S+1 

i=0 ''po •'PO 
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It is enough to show that the system of equation 
s 

(56) Y,Mji(f>,= ^j, j = l,2,...,^ + l 

i=0 

where 

Mji = j ' x'Q{xY^W{x) dx and G 

has a unique solution </)q, . . . , (^5. 

The hnear system (!56|) has N{S + 1) equations and N{S + 1) un- 
knowns, and by Lemma[9]the only solution of the homogeneous system 
is the zero solution. Therefore the system (15^ has a unique solution, 
completing the proof of LemmaO under the assumptions of §4.4[ 

4.5. Proof of Lemma[5] in the general case: analytic and alge- 
braic methods complementing each other. Lemma[T5] shows that 
there is an affine transformation of the dependent variable mapping 
(ITSll into a similar system which has the eigenvalues of the matrices Bj 
increased by one: 

Lemma 13. Consider a nonhomogeneous Fuchsian equation l[T^) . There 
exists a polynomial jq^x) of degree at most 5 + 1 so that the substitution 

(57) y = yo(x) + Q{x)y 
takes pR) into the system 

(58) f{x) + B{x) y(x) = ^ g(x) 
with 

(59) B{x) = B{x) + /, therefore Bj = B^ + I 
and g(x) analytic on D. 



Proof of Lemmau. 
Decompose g as 

(60) g(x) = gp(x) + Q{x)g{x) 

where gp is a polynomial degree at most 5*+ 1 and g is analytic on D. 
Substituting ([^ in ([18]) we obtain ([58]) with 

g(x) = g(x) - yo(x) + gp(x) - B{x) yo(x) 
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The polynomial yo(a;) is now determined so that the difference Sp{^)~ 
B{x) yo{x) has no poles, therefore it is a polynomial in x, degree at most 
S. Indeed, if gj are the residues of gp/Q at x = pji 

/s+i 

(61) g,(x) = Q{x) 



X — Vj 

0=0 



then take 




(62) y„w = ew| >;^-J_^s-., 

giving fl57|) and completing the proof of Lemma[T3l □ 

Remark 1. Note that yo depends linearly on gp, and it does not depend 
on g. 



By repeated use of LemmafTSl we can link solutions of systems ( |T8l) 
with matrices Bj with non-positive real parts to solutions of systems 
with positive real parts. For such systems the existence of corrections 
was proved in §4.4[ What is still needed is to show that these correc- 
tions can be carried back to the initial system. The following Lemma 
shows that this is indeed possible. 

Lemma 14. Consider an equation fT^) and its transformation [58\} 
provided by LemmalTSi 

Let (p be a polynomial of degree at most S and consider the following 
modification of ^5^: 

(63) y'(x) + B{x)y{x) = \g{x) - 0(x) 

Then for any (p there exists 4> a polynomial of degree at most S, so 
that equation l[2^) transforms to using an affine substitution 



(64) y = y^x) + Q{x)y 
with y(f,{x) a polynomial of degree at most S + 1. 

Proof of Lemmall4\ 

It is easy to see (by Remark[T], or by direct calculation) that the 
substitution flMl) transforms fl22|) to fl63|) if y</) = yo — yi with yo(a;) as 
in the proof of LemmafTSl (see ( !60l) . ( !6T1) . ( !62|) ). and yi(x) a polynomial 
solution, degree at most 5 + 1, of 

(65) y'i + i?(,T)yi = l[g^-0 
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The existence of cf) and of the polynomial solution of (l65l) were proved 
in M □ 

Reduction of the general case to the case of §4.41 

Let —n (with n G N) be a lower bound for the real parts of the 
eigenvalues of all the matrices Bq, . . . , Bs+i- 

If n = this means that §4.41 applies to the system. 

If = 1 then using LemmafTSl we obtain a system to which §4.41 
applies, a correction (p is found, and using LemmafT^ a correction (p is 
obtained for the original system. 

If n > 2 then one uses LemmafT^ repeatedly, n times, then §4.41 
applies, a correction is found, and using LemmafTH repeatedly, n times, 
a correction is found for the original system. 

The proof of Lemma[5] is now complete. 



5. Appendices 

Sections §5.H - §5.2I summarize a few known facts regarding solutions 
of linear, first order differential equations, near an isolated singularity 
of the equation. For proofs see p], followed here except for notation. 

5.1. Singularities of linear differential equation. Consider a lin- 
ear system 

(66) ^ = M(x)y, yeC^ 

ax 

If the matrix M{x) is analytic at xq then this point is called a regular 
point of the equation (166!) . It is well known that any first order systems 
(linear or nonlinear) are analytically equivalent to each other near a 
regular point. 

If the matrix M{x) is not analytic at xq then xq is called a singular 
point of the equation. 

If the singular point Xq of M{x) is isolatec^ then every fundamental 
matrix Y{x) of (IMj) has the form 

(67) Y{x) = (!>{x){x-xof 

where $(x) has an isolated singularity at Xq and P is a constant matrix. 

Classification based on the regularity of solutions. If ^{x) has at most 
a pole singularity at Xq, then the point Xq is called a regular singularity. 
Otherwise it is called an irregular singularity. 



^By an isolated singularity it is meant here that M{x) is analytic in a punctured 
disk < |a; — a;ol < 5. In particular M{x) is single- valued at xq- 
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Note that if xq is a regular singularity then, from ( 167|) . the funda- 
mental matrix Y{x) has a convergent expansion as a sum of series in 
integer powers of x — Xq possibly multiplied by noninteger powers of 
X — Xq and integer powers of ln{x — Xq) . 

Classification based on the regularity of the equation. Assume the 
isolated singularity xq of M{x) is a pole: M{x) = (x — xo)~'^~^L{x) 
with L{x) is analytic at x = xq and r > 0. The number r is called the 
rank of the singularity. 

If r = the point Xq is called a singular point of the first kind, or 
a Fuchsian point. It turns out that a Fuchsian point is necessarily a 
regular singuarity. 

More details about the structure of solutions at a Fuchsian point are 
discussed in the following §5.21 



5.2. Systems near a Fuchsian point. Consider a linear system 

(68) ^ = ^—L{x)y, yeC' 

ax X — Xq 

with L{x) a matrix analytic at Xq. 

The structure of the matrix P in fl67p is particularly simple if the 
eigenvalues of L{xq) are linearly nonresonant, in the sense that 

(69) any two different eigenvalues do not differ by an integer 

In this case we can take P = L{xo) and choose a fundamental matrix 
in the form 

(70) Y{x) = ^{x) {x - xo)^^'"'^ with ^{xo) = I 

and $(x) analytic at xq. 

In linearly resonant cases (when there are distinct eigenvalues of 
L(xo) which differ by an integer) the matrix P may differ from L{xq). 
However, the eigenvalues of P (not counting their multiplicity) form a 
subset of the eigenvalues of L(a;o)0. 

Another representation of the general solution in resonant cases is 
again (1701) . only that this time $(a;) is a convergent power series in 
both X — Xq and ln(x — xq). 

Note that in (!35|) there is a different order of multiplication in the 
differential equation, and the Fuchsian point in (!36l) is xq = C)0. 



'''Any resonant string of eigenvalues of L{xa) is replaced in P by the lowest 
eigenvalue in the string with higher multiplicity, see in [3] the Lemma of Chap. 4, 
Sec.4. 
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5.3. Nonlinear perturbation of first order systems with a non- 
resonant Puchsian point. Consider a nonlinear equation 

(71) ^ = ^—L{x)u+^—i{x,u), ueC"" 
ax X — xo X — xq 

with {{x, u) analytic at x = xq, with a zero of order two at u = and 
L(xo) linearly nonresonant. 

Let $(a;) be as in fITOl) . The substitution u = $(x)u in fl7T]) yields 

(72) ^ = -^—L{xo)u + ^^$(x)-if (x, ^x)u) 
ax X — Xo X — Xo 

which is an equation with the same regularity as (I7T1) for small x — xq 
and u and with a constant matrix in the linear part. 

5.4. Proof of Lemma[6l This result appears in pj Ch.5, §22, where 
it is proved for M diagonal, and it is stated for M in Jordan normal 
form. A proof is included here for completeness. 

I. Note that (after a change of coordinates in C^) we can assume 
that M is in normal form. 

Indeed, upon a change of coordinates in C^, taking w to Sw, a C^- 
valued polynomial q(w), q G P„ becomes 5'~^q(S'w). Denote by S# 
this induced change of coordinates on Vn- the operator Jm becomes, 
in the new coordinates, S^^JmS#, which equals Js-^ms (as seen after 
a short calculation). 

II. Note that if M invariates the span of {e^; s < k < t} (for some 
1 < s < t < d) then Jm invariates the span of {r^ j^; s < k < t,\m\ = 
n} - see (1251) for notation. 

This follows because (dwrm,fc)^w is a (polynomial) multiple of ek- 

(73) (dwrm,fc)Mw = 0m(w)efc 

III. Assume M is in normal Jordan form. This means that there is 
a partition 1 = to < ti <...< tp = d so that 

(74) Msj = XjQj + (Tjej^i 

with \j equal for all j with ti < j < tj+i, and with cr^- = 0, and all 
other (Tj = 1. 

Using (I71|) in fl75|) we obtain the formula for the polynomial (j): 

d-l 

(75) (j)^{w) = (A ■ m)w™ + J2 (^j+imjw'^-^^+^^+' 

i=i 

IV. The remarks above imply that the matrix of Jm in the canonical 
base (l28ll is block diagonal. Moreover, with an appropriate ordering of 
this base the matrix of Jm is upper diagonal. 
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This ordering is as follows: we have (m', k') ^ (m, k) if k' < k, or 
k' = k and m' < m in lexicographic order. 

As a example, take d = 3. For A; = 1, we have in increasing ordei@: 
(0,0,n), (0,l,n-l), (0,2,n-2),.., (0,n,0), (l,0,n-l), (l,l,n-2), 
. . . (n, 0, 0), after which comes the same sequence with k = 2 etc. 

The fact that the matrix of Jm is upper diagonal is expressed by the 
fact that Jm i"m,fc is a linear combination of r^i^k' with (m', k') ^ (m, k). 

To show this, consider the first block for M: k from 1 to ti. 

For k = 1 the first vector of the canonical basis is r(o....,o,n).i- From 
we have 0(o,...,o,n)(w) = (A ■ m)w™ and from (1711) it follows that 
Mr(o,...,o,n),i = Aiw'^ei hence r(o,...,o,n),i is an eigenvector of Jm corre- 
sponding to the eigenvalue A ■ m — Ai. For other values of m it follows 
from fl73|) . fl7^ . fl75|) that JM^m,! is a linear combination of the vectors 
rm,i (with coefficient A ■ m - Ai) and r^^ej+ej+ui, j = l,...,d - 1. 
We obviously have (m — ej + e^+i, 1) ^ (m, 1). (Note that the index 
m — Cj + Bj+i is not defined if rrij = 0, or if m^+i = n. Since in both 
cases rrij = these indexes do not appear in fl75l) .) 

For k = 2, . . . ,ti the discussion is similar, only, due to fl74l) also 
earlier terms, with k = 1, appear. 

For the other blocks of M the proof is similar. □ 
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